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ANNUAL STATUS REPORT FOR NASA GRANT NAGW-2013 


The project entitled “A High Resolution Liquid Xenon Imaging Telescope for 
0.3-10 MeV Gamma- Ray Astrophysics: Construction and Initial Balloon Flights” 
was initiated with the submission in August 1992 of the proposal CAL-2015 in 
response to NASA Announcement of Opportunity NRA 92-OSSA-8. 

For the second year of the project, covering the period 1 January 1994 through 
31 December 1994 funds in the amount of $250,000 were allocated, under NASA 
Grant NAGW-2013. 

The results achieved with a 10 liter LXe-TPC prototype during the second year 
can be summarized as follows: 

• The charge response has been measured to be stable to within 0.1% for a period 
of time of about 100 hours. 

• The electron lifetime has been measured to be in excess of 1 ms, corresponding 
to attenuation an length of the order of meters. 

• The energy resolution has been measured to be consistent with previous results 
obtained with small volume chambers. At a drift field of 1 kV/cm the FWHM 
of 662 keV gamma-rays is 12% after corrections for electronic noise, shielding 
inefficiency and rise time effects. At 4 kV/cm the resolution improves to 7% 
FWHM. 

• 3-D gamma-ray imaging has been demonstrated with a non-destructive orthog- 
onal wires readout. A spatial resolution of order 1 mm has been achieved. 

• Monte Carlo Simulation Results: expected background count rate at balloon 
altitude; background reduction algorithms; telescope’s response to pointlike 
and diffuse sources. 

The results obtained with the 10 liter LXe-TPC represent a significant progress 
towards the goal of construction and balloon-flight test of a LXe-Coded Aperture 
Telescope for MeV 7-ray astrophysics. 

Most of the research work carried out during this second year of the project is 
summarized in the three articles attached as Appendix 1—3. 




Appendix 1 



DETERMINATION OF THE GALACTIC 26 AL SOURCE DISTRIBUTION 
WITH A NOVEL IMAGING TELESCOPE 


E. Aprile 

Physics Department and Columbia Astrophysics Laboratory 
Columbia University 

538 West 120th Street, New York, NY 10027 
E.L. Chupp 

Physics Department and Institute for the Study of Earth, 
Oceans, and Space 

University of New Hampshire, Durham NH 03824 

A. Bolotnikov 

NASA/Marshall Space Flight Center 
Huntsville, AL 35812 

P.P. DUNPHY 

Physics Department and Institute for the Study of Earth, 
Oceans, and Space 

University of New Hampshire, Durham NH 03824 


Submitted to: The Astrophysical Journal 
1994 November 14 

COLUMBIA UNIVERSITY 

DEPARTMENTS OF 

PHYSICS and ASTRONOMY 

NEW YORK, NEW YORK 10027 


Columbia Astrophysics Laboratory Contribution Number 557 



DETERMINATION OF THE GALACTIC 26 AL SOURCE DISTRIBUTION 
WITH A NOVEL IMAGING TELESCOPE 


E. Aprile 

Physics Department and Columbia Astrophysics Laboratory 
Columbia University 

538 West 120th Street, New York, NY 10027 
E.L. Chupp 

Physics Department and Institute for the Study of Earth, 

Oceans, and Space 

University of New Hampshire, Durham NH 03824 
A. Bolotnikov 

NASA/Marshall Space Flight Center 
Huntsville, AL 35812 

P.P. Dunphy 

Physics Department and Institute for the Study of Earth, 

Oceans, and Space 

University of New Hampshire, Durham NH 03824 
ABSTRACT 

An instrument is described which will provide a direct image of 7 -ray line or continuum sources 
in the energy range 300 keV to 10 MeV. We illustrate the use of this instrument to study the celestial 
distribution of the 26 A1 isotope by observing the 1.809 MeV deexcitation 7 -ray line. The source 
location accuracy is 2' or better. The imaging telescope is a liquid xenon time projection chamber 
coupled with a coded aperture mask (LXe-CAT). This instrument will confirm and extend the 
COMPTEL observations from the Compton Gamma-Ray Observatory (CGRO) with an improved 
capability for identifying the actual Galactic source or sources of 26 Al, which axe currently not 
known with certainty. Sources currently under consideration include red giants on the asymptotic 
giant branch (AGB), novae, Type lb or Type II supernovae, Wolf-Rayet stars, and cosmic-rays 
interacting in molecular clouds. The instrument could also identify a local source of the celestial 
1.809 MeV 7 -ray line, such as a recent nearby supernova. 


Keywords: nucleosynthesis — 7 -ray astronomy — imaging 7 -ray telescope 
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1. Introduction 

Observation of 7-ray lines from cosmic ray sources provides crucial information to enhance 
our understanding of the origin of the elements. Gamma-ray lines result from radioactive nuclei 
produced in stellar and cosmic ray nucleosynthetic processes as well as by direct nuclear excitation 
of stable nuclei. The first celestial radionuclide identified from its 7-ray emission was 26 A1 (1.03 
Myr meanlife) which decays by positron emission or electron capture to 26 Mg giving a line at 
1.809 MeV. Following the independent predictions of Arnett (1977) and Ramaty and Lingenfelter 
(1977), this line was observed by the high energy resolution spectrometer on the HEAO 3 satellite 
(Mahoney et al. 1984), and was subsequently confirmed by the SMM 7-ray spectrometer (Share et 
al. 1985). Since these satellites detected the largest flux from the direction of the Galactic center 
(~ 5 x 10~ 4 photons cm" 2 s" 1 rad" 1 ), it was tacitly assumed that the putative sources of 26 Al 
are at a typical Galactic center distance of ~ 8.5 kpc. This would imply that ~ 3 - 4 M© of 
6 A1 currently resides in the interstellar medium (ISM). Additional observations of the 26 Al 7-ray 
line made with several balloon borne instruments (von Ballmoos, Diehl, and Schonfelder 1987; 
MacCallum et al. 1987; Varendorff and Schonfelder 1991; Teegarden et al. 1991; Malet et al. 1991), 
confirmed its intensity in the direction of the Galactic center. The instruments used for all these 
observations were telescopes with wide fields-of-view (FOV) and poor angular resolution. 

The COMPTEL and OSSE instruments on the Compton Gamma Ray Observatory (CGRO ) 
are continuing detailed imaging observations of the spatial distribution of the 26 A1 7-ray emission 
with greatly improved sensitivity and angular resolution. The ongoing COMPTEL observations 
have been described in a succession of papers by Diehl et al. (1993a, b, c, d; 1994a, b, c). The 
COMPTEL 1.809 MeV images, currently available for the Galactic plane, reveal broadly peaked 
emission roughly compatible with the molecular (CO) gas distribution, but clumpy emission patches 
are evident at other Galactic longitudes (e.g., Diehl et al. 1994b, c). The central emission, while 
apparently nonuniform, appears to extend over a Galactic longitude range --35° < £ < 40° but 
the maximum is not at the Galactic center. The integrated flux from this inner region of the 
Galaxy is 1.9 x 10 4 photons cm -2 s _1 . Emission enhancements have also been reported from 
an extended region (~ 10° across) in the direction of Vela {£ = 263°) and toward the Eta Carina 
galactic arm (£=285°). The total flux reported for the Vela feature is 2.2xl0 -5 photons cm -2 
s . In the case of emission from the Vela region, two plausible sources (Oberlack et al. 1994) 
are the supernova remnant (£ = 263.9°, b= -3.3°) and the Wolf-Rayet binary WR II (£ = 262.8°, 
b = ; 7 - 7 °)- It is clear fro m these observations that our current knowledge of the distribution of 
Al is dramatically different from that based on the earlier ( pre-COMPTEL ) satellite and balloon 
observations. Basically it was then believed that quasi continuous spatial 7-ray emission would 
reflect the smeared locations of explosive nucleosynthetic events over the past million years in the 
central part of the Galaxy. Now there is a distinct possibility that the required amount of 26 Al 
could be well below the few M© needed for sources at the Galactic center distance if it should be 
confirmed that several relatively close sources such as the (unconfirmed) Vela remnant at a distance 

less than 500 pc contribute. In fact, Diehl et al. (1994c) now estimate that the total Galactic mass 
of 26 Al is <1 M©! 

Further advances in our understanding of the origin of the 26 Al requires instruments with 
higher angular resolution and greater flux sensitivity. Ultimately, instruments capable of high 
energy resolution measurements will be needed to study the profile of the intrinsically narrow line 
™toch could reveal the temperature and motion of the emitting regions (Ramaty and Lingenfelter 

In this paper we first review, briefly, the candidate sources for the 26 Al and the current 7-ray 
observations obtained by COMPTEL (§2). We then present the new imaging telescope (§3.1), with 
a special emphasis on the application of a liquid xenon ionization chamber as three-dimensional 
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(3-D) 7 -ray imager and calorimeter (§3.2), followed by a short summary of recent experimental 
results obtained with a 10 liter liquid xenon prototype detector (§3.3). Monte Carlo simulations 
are presented in §4 to demonstrate the imaging response of the instrument both as balloon- and 
satellite-borne instrument. The instrument’s capabilities will be confirmed by observations of the 
Crab Nebula in a planned balloon borne experiment. 


2. Sources of 26 AL Production 


Whatever the site of production of 26 Al, the basic nucleosynthetic process involved is 25 Mg(p, 7 ) 
26 A1 with the seed nuclide 25 Mg arising from a sequence of reactions, in either explosive C-N burning 
or convective He burning giving the net process 14 N— > 25 Mg. In order for the 26 Al to be visible 
through its decay 7 -ray, it must be carried or ejected from its place of production to optically thin 
regions in a time scale short compared to its halflife. In addition, destruction of Al primarily by 
26 Al(n, p) 26 Mg or 26 Al(n,o) 23 Na must be negligible; thus, the conditions for 26 Al formation are 
a proton-rich, neutron- poor environment. Detailed discussions of this problem have been given 
by Woosley (1986), Clayton and Leising (1987), Woosley (1991), Weaver and Woosley (1993). We 
summarize, briefly, the predicted yield of 26 Al to the ISM from several sources that have received 
serious consideration. We also include the recent proposal by Clayton (1994) that heavy cosmic 
ray interactions in molecular clouds are a possible contributor. 

The Galactic distribution of several of the sources to be discussed is poorly known. Initially, 
T rising and Clayton (1985) have modeled the line emission using a distribution which is strongly 
peaked at the Galactic center to explain the HE AO S and SMM observations. Higdon and Fowler 
(1989) have refined this model and calculated the angular distribution of the 1.809 MeV 7 -ray 
for several modeled distributions which Eire broad, extending to more than a radian from the 
Galactic center. All modeled distributions were smooth except one which followed the distribution 
of molecular clouds which shows emission peaks at the Galactic center and at i ~ 35°. 

It is clear that events on a Galactic scale such as novae, AGB objects, and massive stars (see 
below) which have a low intrinsic yield must occur frequently in the million year meanlife of 26 Al 
in order to produce a line intensity corresponding to the COMPTEL result for the Galactic center 
region. In this case the radiation distribution would be smooth. 

It is important to remember, however, that any single nearby (< 1 kpc) source of 26 Al should 
be considered as potentially capable of producing a measurable flux of the 1.809 MeV 7 -ray line. 

2.1. AGB (Red Giants) 

Following the H and He core burning cycles, main sequence stars of initial mass M < 10 M© 
evolve onto the AGB as red giants of radius (100-1000) R© with a hot white dwarf at its center 
(Iben 1989). Norgaard (1980) and Cameron (1984) suggested 26 Al production by what is called “hot 
bottom burning” (HBB) of H below the convective stellar envelope (Paulus and Forestini 1991). 
For stars with M~5-7 M© these estimates gave 2-3 xl0 ~ 5 M© of 26 Al per AGB star (Prantzos 
1991) which is the maximum that could be injected into the ISM. Using the latest nuclear reaction 
cross-sections, recent calculations for AGB stars with M = 3 M© and metallicity Z = 0.03, (Paulus 
and Forestini 1991; Forestini, Paulus, and Arnould 1991), show that ~ 10 -8 M© of 26 Al per AGB 
star could be added to the ISM after 13 thermal pulses. Extending these calculations to stars with 
main sequence masses 1-3 M©, Forestini et al. (1991) estimate an average 26 Al yield of ~0.1 M© 
injected into the ISM every million years. This result depends on several assumptions for the stars 
in this mass range which include a constant thermal interpulse separation and a constant dredge-up 
26 Al mass per pulse, as well as on the initial mass function (IMF) adopted. 
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2 .2. Novae 

Numerous authors have suggested that thermonuclear nova outbursts are likely sites of 26 A1 
production. The current model of a nova outburst assumes that the event is generated by a 
thermonuclear runaway occurring on the surface of a white dwarf in hydrogen rich gas accreted 
from a companion, i.e., a cataclysmic binary system. Recent calculations by Starrfield et al. (1993) 
for accretion onto the surface of a (O-Ne-Mg) low mass (~1 M©) white dwarf estimate a 26 Al yield 
of a few times 10 7 M 0 per nova, although Weiss and Truran (1990) predicted a yield more than 
an order of magnitude lower. It should be noted that it would require but a single local nova at a 
distance of 15 pc to produce the flux level observed by COMPTEL from the Vela region. 

2.3. Supemovae 

The supernovas source of 26 A1 is expected to be from types SNe lb and SNe II. These events 
are believed to result from the gravitational collapse of massive stars (M> 8 Mq) at the end of 
the successive hydrostatic burning stages and may occur in either isolated or binary systems. 
The recent analysis of Thielemann et al. (1993) shows that O, Ne and Mg (the 26 Al seed nuclei) 
originate mostly in hydrostatic burning shells, with the amount ejected increasing rapidly with 
progenitor mass. Woosley (1991), considering in addition explosive neon burning and various 
neutrino processes (Woosley et al. 1990), shows that the 26 Al yield is generally increased over the 
presupernova value after the explosion for progenitor masses of 15-35 M e (see also Weaver and 
Woosley 1993). Prantzos (1993) has recently reviewed calculations of 26 Al production in massive 
stars including that from Wolf-Rayet stars of mass M>40M© (see §2.4 below). If Wolf-Rayet stars 
are the progenitors of Type lb supemovae, then it is expected that still more 26 Al will be produced 
by explosive nucleosynthesis (see, for example, The et al. 1990). Among several other factors, 
the Al yield depends on the presupernova metallicity which is a principal unknown in theoretical 
estimates of 26 Al production, since it depends on the Galactic locale, decreasing with Galactocentric 
radius (cf. Prantzos 1991, 1993). In Wolf-Rayet stars 26 Al may increase with metallicity (Prantzos 
and Casse 1986), but in SNe II it may actually decrease (Woosley and Weaver 1980). The extreme 
values for the 26 Al yield range from 2xl0- 6 M o to 4x10 ~ 4 M Q for progenitor masses of 12 M 0 to 
35 A/©, respectively. 


2.4. Wolf-Rayet Stars 

Another class of massive stars (M>30 M©) are the Wolf-Rayet (WR) stars which show broad 
bright emission lines against a faint continuum. These stars are also characterized by low hydrogen 
and enhanced helium abundances and in some cases overabundance of N, C and O for the WN WC 
and WO subclasses, respectively (Weiss and TVuran 1992). These stars, near the main sequence, 
produce Al through hydrostatic core H burning at temperatures above T c > 35 x 10 6 K (Dearborn 
and Blake 1985). In general, strong radiatively driven winds from the convective contracting core 
blow off a stellar envelope containing 26 Al. If the mass loss rate and extent of the convective core 
are sufficiently great then 26 Al can be injected into the ISM. 

Recent calculations (Prantzos 1991), using the latest measured reaction rate for 25 Mg (p, 7) 
Al, give 26 Al yields of 10 -5 M 0 to ~ 10 -4 M 0 for stellar massses of 40 M 0 to 100 M 0 . An 
average value of ~ 4.5 x 10“ 5 M 0 is found for M> 40 M 0 using a Salpeter IMF with index ct =- 2 5 
These results are used along with an estimate of 26 Al from SNe II to estimate the total contribution 
from massive stars. An analysis by Prantzos (1993) concludes that WR stars can produce at most 
~20% of the Galactic 26 Al required, while SNe II and SNe lb could produce half of it. (See also , 
Meynet and Arnould 1993.) Since these conclusions depend strongly on our knowledge of the WR 
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radial distribution (as well metallicity and the IMF), uncertainties remain substantial. Azzopardi 
(1992) points out that the number of identified Population I WR stars is only 161 compared to 
the estimated number of Galactic 1200 objects, based on surveys of these distinctive stars in other 
galaxies such as the LMC, SMC, and M31. 

2.5. A Nearby Supernova 

As seen from the above discussions, there are many possible sources for 26 Al, but uncertain- 
ties in the absolute yield predictions from nucleosynthesis calculations and the unknown spatial 
distributions of probable sources will make it difficult to interpret the 26 Al sky maps observed by 
COMPTEL in terms of sources distributed on a Galactic scale. In this regard it is important to 
recognize that local source (s) would reduce dramatically the required amount of radioactivity! 

Some possibilities include: A nearby nova or supernova, nearby massive stars in the Wolf- 
Rayet phase, or, less likely, lower mass AGB stars. The first possibility is made likely by the X-ray 
evidence for a local interstellar bubble, presumably the result of a recent ( < 10 5 years) nearby 
supernova (Cox and Anderson 1982; Cox and Reynolds 1987). The production requirements for 
26 Al could be reduced by a factor of 10 4 for a local single source. 

Dearborn and Blake (1985) have considered the possibility that the massive stellar winds from, 
e.g. Wolf-Rayet stars, could have enriched the local ISM with 26 Al. Morfill and Hartquist (1985) 
estimated the production of 26 Al in a “local” supernova superbubble which suggests again that 
~ 3 — 4 M 0 of 26 Al may not be necessary to explain the observed flux of the 1.809 MeV line. More 
recently Blake and Dearborn (1989) have proposed a specific local source of 26 Al, Wolf-Rayets 
and supernovae production in OB stars in the Sco-Cen association. This model predicts a specific 
distribution of 7-rays in the sky and in particular a 7-ray signal at high Galactic latitude. 

2.6. Cosmic Ray/Molecular Cloud Interactions 

The 7-ray lines recently observed in the direction of the Orion molecular clouds (Bloemen et al. 
1994) are presumed to arise from the interactions of heavy cosmic-rays ( Z > 2) with the hydrogen 
and helium of the ISM. Such an explanation requires that the heavy cosmic-ray flux be enhanced by 
a factor of 30 over that in the solar neighborhood. Using this factor, Clayton (1994) has calculated 
the 26 Al abundance in the cloud from 26 Mg and 28 Si for (H,n) and (H, ppn) reactions, respectively. 
He also assumed that the abundances of Mg and Si relative to C - O in the Orion cloud correspond 
to those in the cosmic rays detected at the Earth and estimated that the production rate of 26 Al 
in the Orion complex could be 5.4 x 10 51 atoms Myr^ 1 . For solar material and a cloud mass of 
10 5 M©, the 27 Al content is estimated as 2.6 xlO 56 atoms, so the 26 A1/ 27 A1 ratio in the cloud 
could be ~ 2.1 xlO -5 . Previously Clayton had estimated that a smaller 26 A1/ 27 A1 ratio of ~ 10" 6 
was needed to obtain 2 M© of 26 Al in the entire ISM, approximately half of which is in molecular 
clouds (see Clayton 1994 and references therein). In principle then, cosmic ray spallation reactions 
in molecular clouds may also be considered as a contributor to the observed 26 Al. 


3. Requirements for Future Experiments 

The instruments that have been used to date for determining the distribution of the 1.809 MeV 
7-rays do not directly provide a celestial image of the radiation. In a double-scatter type Compton 
telescope, such as COMPTEL , the incident direction of a 7-ray can only be determined to within 
a cone of half angle equal to the Compton scattering angle and with a resolution of several degrees 
for the conical shell thickness. The limitation on event topology and range of scattering angles 
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imposed by the two independent detectors of such a system results in a small effective area, and 
determining the proper background is a major problem. For example, one method of determining 
the background for the COMPTEL Galactic 26 A1 line maps was based on the average of high 
latitude observations, which would be compromised by any local soure of celestial 26 A1 (§2.5)), 
but other methods were also used. Complex Maximum Entropy or Maximum Likelihood methods 
convolved with the instrument response were finally necessary in order to extract the image from 
the background events. The line intensities in adjacent sky bins are then not truly independent as 
explained by Diehl et al. (1993a). In the case of the nonimaging, balloon-borne, S MM and HEAO 
3 spectrometers, a collimator defines the direction of a 7-ray signal above background to within an 
angle of typically 20°. Therefore, it is clear that a major advance in 7-ray imaging techniques is 
necessary to determine directly the celestial distribution of 26 Al. 

In order to advance our knowledge of this important astrophysical problem we now describe 
a new imaging coded aperture 7-ray telescope which uses a liquid xenon detector for 7-ray spec- 
troscopy and imaging. As will be demonstrated, the xenon detector can provide an ideal com- 
bination of high detection efficiency over a large area, good energy resolution, excellent spatial 
resolution, and event visualization capability which permits excellent background rejection. When 
used with a coded mask, the combination LXe-CAT gives a 7-ray telescope with excellent celestial 
angular resolution. The simultaneous image of the 7-ray sky and the background field with no as- 
sumptions about the source distribution can provide the information needed in order to determine 
the true nature of the source or sources of 26 A1. 

3.1. The Liquid Xenon-Coded Aperture Telescope (LXe-CAT) 

The proposed telescope (Aprile et al. 1992a) is shown schematically in Figure 1 along with a 
list of its principal characteristics. It combines a liquid xenon time projection chamber (LXe-TPC) 
as 3-D position sensitive 7-ray detector with a coded aperture mask to achieve precise measurement 
of the energy and angular distribution of 7-ray sources in the 0.3-10 MeV energy region. The coded 
mask is a 43x41 element uniformly redundant array (URA) of the type discussed in Fenimore 
and Cannon (1978). The size of the mask elements is 0.81 x 0.85 cm 2 with 1.2 cm thickness. 
With a detector-mask separation of 100 cm, the angular “pixel” is 0.46°x0.48° over a 19.3° x 19.3° 
FOV. The telescope angular resolution is 30 arcminutes but point sources can be located with an 
accuracy as good as a few arcminutes, based on the excellent signal-to-noise ratio of the liquid 
xenon imaging detector. A larger mask-detector separation and/or smaller mask elements will give 
even more precise source location. The LXe detector is surrounded by an active anticoincidence 
shield to reduce the cosmic and atmospheric background 7-ray flux at balloon altitude. The active 
area of the LXe-TPC is 1200 cm 2 . With an active liquid xenon layer of 10 cm ( Z — 54, density = 
3.06 g cm -3 ), the full energy peak efficiency is about 65% at 1 MeV. 

The LXe-TPC works on the principle that the free ionization electrons liberated in the liquid 
by a charged particle can drift under a uniform electric field, inducing a signal on sensing electrodes 
(see Figure 2). For a 7-ray, it is the electron resulting from a photoelectric absorption or a Compton 
scattering, or the electron-positron from a pair production, which will ionize and excite the xenon 
atoms yielding a large number of electron-ion pairs and scintillation photons (Doke et al. 1990, 
Kubota et al. 1977). The scintillation light signal is fast (< 10 ns), and thus an ideal marker 
of the time origin of a 7-ray interaction. From the measured electron drift time, referred to this 
time origin, and the known drift velocity, the coordinate of the interaction point, along the drift 
direction, is inferred. The other two coordinates, as well as the 7-ray energy, are inferred from the 
charge signals induced on the sensing electrodes. Thus the LXe-TPC measures both the energy and 
the spatial distribution of each ionizing event occurring within the sensitive volume. The factors 
which determine the ultimate energy and spatial resolution of a liquid xenon chamber have been 
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studied. Experimental results show that a submillimeter spatial resolution and a 6% FWHM energy 
resolution at 1 MeV can be realized in a practical detector (Aprile et al. 1991, 1992c). Improved 
energy resolution has been measured with liquid xenon doped with appropriate photosensitive 
materials (Ichinose et al. 1992; Aprile 1994a). 

The intrinsic 3-D imaging capability of the LXe-TPC allows a direct visualization of the com- 
plex histories of 7-ray events with multiple interactions, initiated by either Compton or pair produc- 
tion, thus providing a powerful method for background identification and rejection. Using Compton 
kinematics, the original direction of a 7-ray undergoing multiple Compton scatterings in the liq- 
uid xenon can in fact be reconstructed from the measured energy and spatial coordinates of each 
scattering point. Events which kinematically could have not have come from the telescope’s FOV 
are rejected as background. The signal-to-noise ratio can be further improved through the iden- 
tification and rejection of background 7-rays interacting via photoabsorption or pair production, 
based on their location within the sensitive volume and their topology. The background rejection 
capability of the LXe-TPC, by Compton event reconstruction and by fiducial volume selection, 
has been demonstrated with Monte Carlo simulations (Aprile et al. 1993) and, more recently, with 
experimental results obtained with a 10 liter laboratory prototype (Aprile et al. 1994b). 

The response of the LXe-TPC as a Compton polarimeter has also been studied, showing good 
sensitivity to polarization fractions as small as a few percent for strong sources, even in a balloon 
flight (Aprile et al. 1994a). 

3.2. The LXe-TPC as 3-D Position Sensitive Detector 

To image the ionization electrons produced by 7-ray interactions inside the sensitive volume of 
the LXe-TPC, we have proposed the system of sensing electrodes shown in Figure 2. The system 
consists of two orthogonal sense wire planes, separated from the drift region by a screening wire 
grid, and followed by a segmented plate anode. The presence of the grid means that a signal is 
induced on the sense wires only once the drifting electrons have crossed the grid. 

The induction signal is positive while the electrons move toward the wires and negative while 
they drift away from the wires. The positions of the wires hit gives the X — Y coordinates of 
the electron image. The drift time, referred to the zero time of the scintillation light trigger, 
plus the known drift velocity, provides the 2-coordinate information. The total event energy is 
measured from the charge collected on the plate anode. The idealized pulse shapes for a point-like 
charge moving in the drift region are also shown in Figure 2. While the induction signal on the 
wires is triangular in shape, the collection signal on the anode is step-like. For 'y-v&y events with 
multiple Compton interactions, the anode signal has a characteristic multi-step signature. The 
height of each step is proportional to the the energy liberated in the interaction point, and the time 
interval between the steps is proportional to the relative distance along the 2-direction between 
the interaction points. The sum of all the step pulse heights is proportional to the total energy of 
the original 7-ray, if it is totally absorbed. 

The readout structure shown in Figure 2 is ideal for imaging the point-like charge clouds 
produced by MeV 7-rays in liquid xenon with a spatial resolution better than d f \/l2 , where d = 2 
mm is the sense wire spacing. This is due to the strong dependence of the induction signal amplitude 
on the lateral position of the charge cloud with respect to the sense wire (Gatti et al. 1970). The 
center of gravity of the image can then be determined from a weighted average of the signals induced 
on two adjacent wires, with an accuracy which depends on signal-to-noise. Other factors which 
influence the spatial resolution of the LXe-TPC are: 1) the physical extension of the charge cloud 
in the liquid, and 2) its diffusion during the drift. The intrinsic size of the charge cloud produced 
by 1 MeV ionizing radiation in liquid xenon, estimated by Monte Carlo, is on the order of 0.5 mm, 
which is not negligible compared with the detector’s granularity. On the other hand, the spread 
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by diffusion of the center of gravity of a cloud of N electrons is given by ^/2Dt d /N, where D is 
the diffusion coefficient and t d is the drift time. Taking D — 65 cm 2 s _1 (Shibamura et al. 1979), 
t d = 40/rsec as maximum drift time in the proposed LXe-TPC, and N = 6.4 x 10 4 for 1 MeV 
ionizing radiation, the resulting diffusion spread is only 3 /im. 


3.3. Test Results with a 10 Liter LXe-TPC Prototype Detector 

To demonstrate the operation of a large liquid xenon detector, particularly spectroscopic and 
3-D imaging capabilities, we have built and are currently testing a 10 liter LXe-TPC prototype, 
implemented with the electrode system discussed above (see Aprile et al. 1994b, Chen 1994). The 
prototype detector has half the drift gap, one third the sensitive area, and a volume similar to the 
sensitive volume of the proposed flight instrument. A spacing of 4 mm for theX-Y wire planes was 
chosen to minimize the total number of electronics channels. Experiments have been carried out to 
test the cryogenics operation, the level of liquid xenon purity (which affects charge collection and 
thus energy resolution) and the imaging response to various 7-ray sources. We have verified that 
10 liters of liquid xenon can be purified and maintained with a contamination level < 1 ppb over a 
period of several days. We have reproduced the energy resolution results previously obtained with 
smaller and simpler liquid xenon chambers. For the first time we have simultaneously measured the 
energy and the image of a collimated beam of 7-rays. Figure 3 shows the reconstructed image of 
662 keV 7-rays from a 137 Cs source collimated with a 2 mm diameter hole. The data are consistent 
with a detector RMS spatial resolution better than 1 mm. 

Event reconstruction based on Compton kinematics has also been demonstrated. A parallel 
vertical beam of 7-rays from a point source was simulated by placing a radioactive 7-ray source 
40 cm above the detector’s sensitive area. Events with a single Compton scattering followed by a 
photoabsorption were selected for this analysis. As an example, Figure 4 shows the on-line display 
of a 1.274 MeV 7- ray event from a 22 Na source-beam. The magnified view of the signal on the 
anode clearly shows a two-step event indicating that the 7-ray history was a Compton interaction 
followed by photoabsorption. The sum of the two steps pulse heights corresponds to the total 
energy of 1.274 MeV. From the time and amplitude analysis of the induction and collection si gnals , 
the coordinates and the energy for the two interaction points are inferred as well as their spatial 
separation. The most probable scattering angle is then found from the kinematics of the Compton 
scattering process. The accuracy on the angle determination depends on the spatial and energy 
resolution as well as on the separation between the two consecutive interaction points. Figure 5 
shows the reconstructed event circles for nine 22 Na 7-ray events with a Compton scattering followed 
by photoabsorption. These experimental results directly demonstrate the LXe-TPC capability to 
identify and reject background via Compton event selection. In addition, they verify the feasibility 
of the original proposal (Aprile et al. 1989) to use a LXe-TPC as Compton/pair telescope for high 
energy 7-rays. 


4. Monte Carlo Simulations 

To demonstrate the capability of the LXe-CAT as an imaging 7-ray telescope, we have carried 
out Monte Carlo calculations that simulate the response of the instrument to several source dis- 
tributions. In particular, we simulate the imaging process for a coded aperture imaging system in 
response to photon continuum spectra from point sources and to a 7-ray line from localized sources. 
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4.1. Simulation of Coded Aperture Imaging 

Coded aperture imaging is related to the “multiple pinhole” technique (Dicke 1968). A planar 
absorber with a pattern of openings is placed between a 7 -ray source and a position-sensitive 
detector. The absorber (or aperture) modulates the incident 7 -ray flux, producing a response 
pattern on the position-sensitive detector. In general, the pattern is not a recognizable image 
of the source, so the pattern must be decoded or deconvolved to produce an image. Ideally, the 
imaging response to a point source should be a delta function with flat sidelobes. It has been shown 
(Fenimore and Cannon 1978) that a class of aperture patterns called Uniformly Redundant Arrays 
(URAs) has these properties. The source image can be reconstructed from the measured detector 
response pattern using a simple correlation calculation. The response pattern is correlated with 
a decoding function that is closely related to the aperture pattern. The correlation procedure is 
simply 

S' = P*G (1) 

where S' is the reconstructed image, P is the detector response, G is the decoding function, and * 
is the correlation operator. In terms of discrete arrays, 


S\i, j) = ^E P(k, l)G(i + k,j + 1) (2) 

k = 1 /=1 

where r and s are the dimensions of the (rectangular) aperture pattern. A review of the properties 
of coded aperture imaging can be found in Caroli et al. (1987). 

The LXe detector is simulated as an active volume of liquid xenon (density — 3.06 g cm -3 ) 
that is 35 cmx35 cm and 10 cm thick. The calculations also simulate inactive xenon around the 
sensitive volume in the detector container, and a 5 cm thick Csl shield surrounding the LXe detector 
and extending 20 cm above the top face. The coded aperture mask is modeled as a URA pattern 
of blocks of tungsten alloy (90% W and 10% Cu, density = 17.0 g cm“ 3 ). The exact pattern that 
is used depends on the desired angular resolution^ 

For our simulations, we choose a mask design and a mask/detector separation that will give 
sub-degree angular resolution over a reasonably wide field-of-view (FOV). In this case, the basic 
mask pattern is a 43 x 41 element URA at a distance of 100 cm from the LXe detector face. 

Monte Carlo histories are generated for 7 -rays incident on the telescope at energies between 
0.3 and 10 MeV. Photoelectric, Compton, and pair events are identified and recorded. Secondary 
electrons, resulting from 7 -ray interactions, are assumed to lose all their energy at the point of 
creation if their kinetic energy is less than 2 MeV. The program also includes the uncertainties in 
the interaction location and in the energy loss at the interaction point due to the detector’s finite 
spatial and energy resolution. 

A point source of 7 -rays is simulated by a parallel beam of photons that is uniformly distributed 
over the active area of the detector. The program checks for interactions in the mask and the active 
xenon volume, as well as in other sections of the instrument. The program records the location of 
interactions and follows the histories of the secondary 7 -rays and electrons. An inventory is also 
kept of the energy depositions within the active xenon volume. For events with multiple interaction 
points, the algorithm based on Compton kinematics is used to determine which was most likely to 
be the point of the first interaction (e.g., Aprile et al. 1993), since this information is required for 
the imaging process. 

We simulate the background in the shielded LXe detector with a random distribution of counts 
throughout the detector volume. The number of background counts is determined by the detector 
efficiency, the background flux expected in the operating environment, and the observation time. 
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The following components of the continuum instrumental background were considered to estimate 
the total background count rate in the LXe-CAT at balloon altitude: 1) aperture flux; 2) shield 
leakage; 3) elastic neutron scattering; 4) spallation products; and 5) background from the coded 
aperture mask. Cosmic diffuse and atmospheric 7-rays are the major sources that contribute 
to components 1) and 2) at balloon altitudes. The parameterized form of the atmospheric 7- 
ray spectrum given by Costa et al. (1984) and the cosmic diffuse spectrum given by Schonfelder, 
Graser and Daugherty (1977) were used. The neutron spectrum inside the LXe detector was taken 
from Gehrels (1985). For the spallation spectrum, we used the spectrum from Dean, Lei, and 
Knight (1991) scaled to the LXe density and mass number. The background from the tungsten 
mask (due mainly to secondary photons from the interaction of atmospheric 7-rays in the mask) 
has been estimated using the empirical formula of Owens et al. (1985, see also Dunphy et al. 

1989) . Background in the LXe-CAT instrument in a low Earth orbit was also estimated using the 
techniques described in Gehrels (1992). 

Figure 6 shows the calculated total background spectrum and its individual components at 
3 g cm over Palestine, Texas. To reduce the background, we apply the event reconstruction 
algorithm based on Compton kinematics discussed in Aprile et al. (1993). At 2 MeV the estimated 
count rate is reduced by a factor of ~ 2 when events which kinematically could not have entered 
the aperture defined by the shield ( ~ 50°) are rejected. By applying a more stringent software 
collimation of only 3°, for the case of a very narrow FOV, the background is reduced by almost 
an order of magnitude. For a detailed description of the background calculations and the effective 
reduction of background events in the MeV region we refer to Aprile et al. 1994c. 

In general, the background in the detector will not be distributed uniformly throughout the 
sensitive volume. If not accounted for, this will cause a false structure to appear in the reconstructed 
image. In practice, this false structure can be removed by making a background observation under 
similar conditions for a time equal to the source observation time. The background pattern in the 
position sensitive volume can then be used to correct the source observation, cancelling out the 
structure due to the background non-uniformity. Because of the background subtraction, however, 
the effective background (i.e., its variance) is increased by a factor of 2. T his approach has been used 
successfully in balloon flights of a coded aperture 7-ray detector (Dunphy et al. 1989; Bhattacharya 

1990) . 

The source plus background simulation produces a distribution of event locations in the detec- 
tor. The locations are binned to match the dimensions of the mask elements (i.e., anrxs array), 
producing the array P(k, l ) in equation (2). The correlation of this array with the decoding array, 
G, produces an image of the source, S'. 

4.2. Simulation Results 
4.2.1. The Crab Nebula Region 

To illustrate the response of the LXe-CAT to continuum spectra of 7-rays from discrete sources, 
we simulate the Crab Nebula. The differential photon spectrum is a power law of the form (Pen- 
ningsfeld, Graser and Schonfelder 1979) 

dN/dE = 3 x 10- 3 E^ cm-V^eV- 1 . (3) 

First, we simulate images from a typical high-altitude balloon observation. The source observation 
time is taken to be 10 hours, with an equal length of time taken for measurement of background 
non-uniformity. The background shown in Figure 6, reduced by a factor of ~ 2 after Compton 
event reconstruction, is applied. This corresponds to the background rejection expected for a 19° 
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FOV. Absorption of source photons by the residual atmosphere at a depth of 3.0 gem 2 is included 
in the calculations. The effective increase in background (a factor of 2) required by the evaluation 
of background non-uniformity is also taken into account. Images resulting from the simulation in 
two energy ranges are shown in Figure 7. The source is present at a significance level of 16<x in the 
range 0.5-1 MeV (Figure 7a) and Ilex in the range 1-10 MeV (Figure 7b). Since a point source 
can be located with a precision of (McConnell et al. 1987) 

&loc = Gpixell^o 

and the characteristic size of a pixel is about 30', the uncertainty in the Crab Nebula location for 
this observation would be about 2'. 

To illustrate the capability of the same detector to resolve separate point sources in the same 
FOV, we simulate a satellite-based observation of the Crab Nebula and the quasar PKS0528+134. 
For the 7-ray spectrum from the quasar, we use the most intense spectrum > 3 MeV consistent 
with the COMPTEL observations reported by Collmar et al. (1993). The live-time is taken to be 2 
weeks and the background for a low Earth orbit is used. The image in the energy range 3-20 MeV 
is shown in Figure 8. PKS 0528+134 is seen in the image at a significance level of 6<x. 

4.2.2. The Galactic Center and Vela Regions 


As discussed in §1, the measurement of the 26 Al distribution in the Galaxy via the 1.809 MeV 
7-ray line is an extremely important objective for 7-ray astronomy. For the 1.809 MeV source 
distribution, we use the maps of this radiation observed by COMPTEL (Diehl et al. 1994c). These 
maps show a highly structured emission over a wide longitude range with “hot spots” scattered 
along the Galactic plane, especially near the Galactic center and in the direction of the Vela and 
Cygnus regions. 

For the Galactic center “hot spot” we use the flux of 3.6 xl(T 4 cm s reported by Diehl 
et al. (1994b) and assume a “point source” distribution (i.e., less than about 0.5° in extent, the 
intrinsic angular resolution of the simulated detector). We assume a satellite-based observation 
time of 1344 hr. This corresponds to the COMPTEL observation live-time (~ 8 weeks). The 
simulated image for these conditions is shown in Figure 9. The significance of the galactic center 
“hot spot” is about 13<x. This compares favorably with the 6.5<x significance of this source for 
COMPTEL (Diehl et al. 1993). The Vela “hot spot,” with a flux of 2.2 xlO -5 cm -2 s _1 , would 
have a significance of 9<x in the LXe-CAT for the same observing conditions. 

We note that if the source is extended in size such that the flux is divided equally among n 
image pixels, the source significance in each pixel will be reduced by a a factor of n. An improved 
response to extended emission can be achieved by a modified LXe-CAT configuration or by removing 
the coded aperture and using the LXe-TPC as a Compton telescope. 


5. Conclusions 

We have described the design of a unique 7-ray imaging telescope operating in the energy 
range 0.3—10 MeV. The telescope is based on a liquid xenon time projection chamber (LXe-TPC). 
Location and energy analysis of events with multiple Compton scatterings in the liquid permits to 
use Compton kinematics to reconstruct the possible directions of the incident photon. This will 
allow significant reduction of background. Use of fiducial volume cuts will also allow background 
reduction for single scattering events and internal /3-decays. 

Another major benefit of using a TPC is that the heavy external shielding usually added to 
attenuate background can be eliminated or greatly reduced because of the background reduction 
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techniques allowed by the intrinsic 3-D imaging. By combining the LXe-TPC with a coded mask 
we obtain an imaging telescope (LXe-CAT) with superior celestial angular resolution and source 
ocalization accuracy over other instruments. We have chosen to demonstrate the capabilities of 
this new instrument by performing simulations for a satellite observation of localized sources of 
1.809 MeV A1 identified by COMPTEL. Since it is essential to verify the properties of a new 
telescope design with a balloon test flight, we have also performed a simulation of a balloon-based 
observation of the Crab Nebula 7-ray spectrum. 

Considerable work has been done at the Columbia University Astrophysics Laboratory to assess 
the feasibility of a liquid xenon chamber for 7-ray imaging and spectroscopy. A 10 liter LXe-TPC 
prototype has been developed and tested in the laboratory. A first balloon flight of a LXe-CAT to 
observe the Crab Nebula and verify the instrument’s imaging capability and the background in a 
space environment could be carried out within two years. 
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Figure Captions 


Figure 1. Schematic view of the Liquid Xenon-Coded Aperture Telescope (LXe-CAT). 


Figure 2. Schematic view of the electrodes system for the three dimensional Liquid Xenon-Time 
Projection Chamber (LXe-TPC). 


Figure 3. Two-dimensional spatial distribution of 662 keV 7-rays from a collimated 137 Cs source 
as detected by the 10 liter LXe-TPC prototype. The collimator diameter is 2 mm. 


Figure 4. The on-line display of a Compton scattering event recorded in the 10 liter LXe-TPC 
prototype, showing the digitized induction signals and anode signal. The anode signal has been 
magnified to show the double interaction signature of the event. 


Figure 5. Reconstructed event circles are shown for nine Compton scattering events from a point 
source of 22 Na 1.274 MeV 7-rays. The celestial sphere is subdivided into 5° x 5° bins and the source 
location is at 0° x 0°. 


Figure 6. The estimated differential background count rate spectrum and its components are shown 
for the LXe-CAT at balloon altitude (3 g cm" 2 ) over Palestine, Texas. The LXe-CAT was shielded 
with 5 cm of Csl and no “software” background rejection was applied. 


Figure 7. A simulated image of the Crab Nebula is shown for a 10 hr balloon observation, a) image 
in 0.5-1 MeV energy range (16a significance); b) image in 1-10 MeV range (11a significance). 


Figure 8. A simulated image of the Crab Nebula and quasar PKS 0528+134 in the 3-20 MeV 
energy range is shown for a 336 hr satellite observation. PKS 0528+134 is significant at a level of 
6a. 


Figure 9. A simulated image of the 26 A1 “hot spot” near the Galactic center at the energy 1.809 
MeV is shown for a 1334 hr satellite observation. The source is significant at a level of 13a. 
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Abstract. We describe a unique 7 -ray imaging telescope operating in the energy range 
0.3-10 MeV. The basic element of the telescope is a liquid xenon time projection cham- 
ber (LXe-TPC) as 7 -ray spectrometer and 3-D imager. Location and energy analysis of 
multiple Compton scattering events in the chamber will permit a reconstruction of the 
incoming 7 -ray direction, thus allowing a significant reduction of background. A 10 -liter 
LXe-TPC prototype, with 400 cm 2 active area and 5 cm drift gap has been built and 
is currently being tested with 7 -ray sources in the laboratory. Initial results are present- 
*ed. The combination of the LXe-TPC with a coded aperture (LXe-CAT), results in a 
telescope with superior detection efficiency, angular and energy resolution, and excellent 
background rejection capability. Simulation results on the instrument’s ability to image 
the Crab Nebula region in a balloon observation are presented, as well as to image the 
1.809 MeV 26 A1 distribution from the Vela region, in a satellite observation. 

Key words: 7 -ray Astronomy - Imaging - Instrumentation 


1 . Instrument Description 

The proposed telescope (Aprile et al. 1992a) combines a liquid xenon time 
projection chamber (LXe-TPC) as calorimeter and position sensitive 7 -ray 
detector with a coded aperture mask to achieve precise measurement of the 
energy and angular distribution of 7 -ray sources in the 0.3-10 MeV energy 
region (Fig. 1 ). The coded mask is a 43x41 element URA, with elements 
0.81x0.85x1.2 cm 3 in size. With a detector-mask separation of 100 cm, the 
angular “pixel” is 0.46° x 0.48° over a 19.3° x 19.3° FOV. The telescope point 
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LXe-CAT Properties 


Energy Range 

0.3 - 10 MeV 

Energy Resolution 

3.9% FWHM at 1 MeV 

3-D Spatial Resolution 

< 1 mm 

Geometrical Area 

**1200 cm 2 

Active Shield 

5 cm Cal (or LXe) 

JFOV (Polly Coded) 

19* x 19* FWHM 

Angular Resolution 

30 arc min 

Point Source 


Location Accuracy 

-2 arcmin(10cr eouroe) 

Kfin Rhx (Line) 

IxUX’phcor 2 * 4 

3oat 1 MeV 

(3 x 10 4 *) 

Min Flux (Continuum) 

3x lC 7 ph cnf 2 » 4 keV 4 

3o«t 1 MeV 

(3x 10 4 *) 


Fig. 1 Schematic view of the Liquid Xenon-Coded Aperture Telescope (LXe-CAT) 

source location accuracy is estimated at 2' for a 10 a source, due to the 
excellent signal-to-noise ratio of the imaging detector. The LXe-TPC has 
an active area of 1200 cm"* and is surrounded by an active anticoincidence 
shield consisting of Csl or LXe scintillator. With a 10 cm deep layer of LXe 
(Z = 54, density = 3.06 g cm -3 ), the full energy peak efficiency is about 
65% at 1 MeV. 

The LXe-TPC works on the principle that the free ionization electrons 
liberated by a 7-ray interaction in the liquid can drift, under a unifor m elec- 
tric field, toward a signal readout structure. The ionization signals induced ' 
on the readout sensing elements provide both the spatial and the total ener- 
gy information for each event. In our design, the information in the X — Y 
plane is obtained from the signals induced on two orthogonal wire planes, 
while the Z information, along the direction of drift, is inferred from the 
drift time, measured with respect to a time zero (see Fig. 2). The fast 
signal (< 10 ns) from the primary scintillation light of LXe, is used for 
time zero mesurement. With the capability of measuring the three spatial 
coordinates and the energy deposited for each 7-ray interaction, the TPC is 
therefore ideal for event reconstruction based on Compton lrinematfea The 
direct outcome of this event visualization is the capability of background 
rejection for both single- and multiple-sites energy deposition events. This, 
as well as the polarization sensitivity of the LXe-TPC as Compton polarime- 
ter, has been demonstrated with Monte Carlo simulation results (Aprile et 
al. 1993b, 1994a). The technical feasibility of such a detector and the factors 
which determine its ultimate energy and spatial resolution have been stud- 
ied at Columbia for the past few years (Aprile et al. 1993a and references 
therein). 
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Fig. 2 Schematic view of electrodes system Rg. 3 On-line display of a Compton scattered y-ny 

for 3-D position sensitivity of LXe-TPC went (1.274MeV), showing the digitized waveforms 

induction signals and anode signal 


2. Test Results with the 10-liter LXe-TPC Prototype Detector 

To demonstrate the operation of a large liquid xenon detector, we have built 
and are currently testing a 10-liter LXe-TPC prototype, implemented with 
the electrode system shown in Fig. 2 (see Aprile et al. 1994b). With a wire 
spacing of 4 mm, the total number of sense wires for X-Y readout is 96. The 
prototype has half the drift gap, one third the sensitive area and a volume 
similar to that of the proposed flight instrument. Experiments have been 
carried out to test the instrument’s spectroscopic and imaging performance. 
The detector’s charge yield was as expected in high purity LXe and stable 
over a maximum period of 100 hr. 

To simulate a parallel beam of 7-rays from a point source, a 60 Co or 22 Na 
source was placed 40 cm above the detector’s sensitive area. Events with a 
single or multiple Compton scatterings, as well as photoabsorption events 
were accumulated. As an example, Figure 3 shows the on-line display of a 
1.274 MeV 7-ray event from the 22 Na source experiment. The magnified 
view of the signal on the anode clearly shows a two-step event indicating that 
the 7-ray history was a Compton interaction followed by photoabsorption. 
The sum of the two steps pulse heights corresponds to 1.274 MeV total 
energy. From the amplitude and time analysis of the induction and collection 
signals, the coordinates and the energy for the two interaction points are 
inferred as well as their spatial separation. The most probable scattering 
angle can then be found from Compton scattering kinematics. These initial 
results directly demonstrate that the LXe-TPC can be very effective for 
background rejection based on Compton event reconstruction. In addition, 
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Fig. 4 2-D imaging of 662 keV ^rays from collimated ,37 C» source 
(2 mm 0 collimation), detected with the 10 liter LXe-TPC 


they verify that the original proposal (Aprile et al. 1989) to use a LXe- 
TPC as a Compton/pair telescope with large efficiency and good angular 
resolution is feasible for high energy 7 -ray astronomy. As an indication of 
the detector’s imaging performance, Fig. 4 shows the reconstructed image 
of a collimated beam (2 mm diameter collimator) of 662 keV 7 -rays from a 
137 Cs source. The data are consistent with a detector RMS spatial resolution 
of about 1 mm. 


3. Simulation Results 

To demonstrate the capability of the LXe-CAT as an imaging 7 -ray tele- 
scope, we have carried out Monte Carlo calculations that simulate the re- 
sponse of the instrument to several source distributions. In particular, we 
simulate the imaging process for a coded aperture imaging system in response 
to photon continuum spectra from point sources and to a 7 -ray line from 
localized sources. The background rates expected in a balloon or satellite 
based observation were calculated. Fig. 5 shows the total background spec- 
trum, together with its individual components, calculated for 3 gem -2 over 
Palestine, TX. A 5 cm Csl thick shield was assumed for the r-a. 1 riiiat.irwi 
To illustrate the response of the LXe-CAT to continuum spectra of 7 -rays 
from discrete sources, we simulate the Crab Nebula for a typical balloon 
observation. The Crab spectrum is taken from Penningsfeld et al. 1979. The 
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Rg. 5 Calculated background spectrum in LXe- 
CAT at 3 g cm' 2 over Palestine, TX. A 5 cm thick 
Csl shield is assumed 


Rg. 6 Simulated image of Crab Nebula (0.5-l.OMeV) 
for 10 hr balloon observation 




Rg. 7 Simulated image of Cmb Nebula and quasar Rf * Simulated image of Vela "hotpot" at l.S09MeV 

PKS 0528+1 34 (3-20MeV) for 336 Ir satellite based (^AI) for as 8-week satellite-based observation 

observation 


source observation time is taken to be 10 hours, with an equal length of time 
taken for measurement of background non-uniformity. The background rate 
of Fig. 5 was reduced by a factor of ~ 2 due to event rejection based on 
Compton reconstruction, but was increased by a factor of 2 as required by 
the evaluation of background non-uniformity. The source is present at a 
significance level of 16cr in the range 0.5-1 MeV (Fig. 6) and 11<7 in the 
range 1-10 MeV. Since a point source can be located with a precision of 
&loc = CTpi xe i/n a (McConnell et al. 1987) and the characteristic size of a 
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pixel is about 30', the uncertainty in the Crab Nebula location for this 
observation would be ~ 2'. 

To illustrate the capability of the detector to resolve separate point 
sources in the same FOV, Fig. 7 shows the image for a satellite-based obser- 
vation (336 hr) of the Crab Nebula and the quasar PKS 0528+134 (Collmar 
et al. 1993). The image of the quasar, in the energy range 3 — 20 MeV, is 
seen at a significance level of 6 a. 

Finally, the measurement of the 26 Al distribution in the Galaxy via the 
1.809 MeV 7 - ray line is an important objective for 7 -ray astronomy and for 
the LXe-CAT instrument (Aprile et al. 1994c). The maps of this radiation 
observed by COMPTEL (Diehl et al. 1994), show a highly structure emission 
over a wide longitude range with “hot spots” near the Galactic center the 
Vela region. Based on these maps, we have obtained the image shown in Fig. 
8 for the Vela “hot spot.” The result is for a satellite observation time of 
1344 hr (similar to the COMPTEL observation time). The Vela “source,” 
with a flux of 2.2 x 10 -5 cm -2 s , would have a significance of 9 a in the 

LXe-CAT. 

This work was supported by a NASA Grant (NAGW-2013) to the Columbia 
Astrophysics Laboratory. 
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ABSTRACT 


The experiment LXe-CAT (Liquid Xenon-Coded Aperture Telescope), which we have proposed for 7-ray as- 
trophysics observations in the 300 keV - 10 MeV energy range, uses a Liquid Xenon Time Projection Chamber 
(LXe-TPC) as a position sensitive 7-ray detector, and a coded aperture mask to provide a telescope with an angular 
resolution of 30' over a Field-of-View (FOV) of 28° x 20°. The point source localization accuracy is about 2' for a 
10<r source. 

To demonstrate the spectroscopy and imaging capabilities of the LXe-TPC as a 7-ray detector we have designed, 
built and are currently testing a 10 liter prototype with a sensitive area of ~ 400 cm 2 and an active liquid depth 
of 5 cm. Studies on stability, spectroscopy and imaging response of this detector to various 7-ray sources will be 
presented. The unique capability of the LXe-TPC to identify and reject background events based on Compton 
kinematics reconstruction has also been demonstrated. 

1. INTRODUCTION 


The study of astrophysical 7-rays in the MeV region has always been challenging from an observational point 
of view, mainly due to the small source flux rate compared with the background rate, the small interaction cross- 
section and the complexity of Compton scattering which dominates in this energy range. Progress in MeV 7-ray 
astronomy relies on new telescopes with superior imaging capability, high energy resolution, wide sky coverage, and 
low background. 

Two important astronomical observations which typify the need for novel imaging instruments are the 1.809 
MeV line emission from radioactive 26 A1 and the 0.511 MeV positron-electron annihilation line emission. The exact 
origin of these emissions is still undetermined. For both cases it is believed that a precise mapping of the spatial 
distribution of the emission will uniquely identify the origin of this 7-ray radiation based on the known distribution 
of the candidate sources 1 . 

Among the proposed novel techniques for imaging 7-ray sources in the MeV energy region, the LXe-TPC is 
recognized as very promising. The properties of LXe make it an excellent radiation detection medium for 7-rays. 
When used in an ionization chamber, operated in the time projection mode, LXe offers an ideal combination of 
high detection efficiency, three-d imensio nal (J*D) event imaging with submillimeter spatial resolution, and very good 
energy resolution response. Like an electronic bubble chamber, a TPC is capable of visualizing the complex histories 
of 7-ray events with multiple interactions, initiated by either Compton scattering or pair production. As a result 
of this imaging, efficient background rejection is also achieved, as well as sensitivity to 7-ray source polarization 
measurements. 

To overcome the limited angular resolution of a LXe-TPC as a Compton telescope for MeV 7-rayB, 2 we have 
proposed the LXe-CAT instrument, shown schematically in Fig. 1, for studies of galactic and extragalactic 7-ray 
sources in the energy range 0.3 - 10 MeV 3 . The telescope consists of a LXe-TPC as position sensitive detector, an 
active anticoincidence shield, to reduce the background rate at balloon altitude, and a coded mask. A list of the 
telescope’s principal characteristics is also included in Fig. 1. The coded mask is a 2 x 2 mosaic of a basic uniformly 
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Fig- 1 Schematic of the LXe-TPC/coded aperture imaging y-ray telescope 


wduniJant array (URA) 4 . It consists of a 43 x 41 dement pattern of 0.91 x 0.68 x 1.2 cm 8 thick tungsten blocks. 
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sensitivity based on a non-destructive signal readout. The results obtained give us the confidence in the performance 
of LXe 7 -ray detectors and their applications in astrophysics. In order to demonstrate the direct imaging capability 
of a large volume LXe-TPC, we have designed and built a 10 liter prototype, which has about one third of the 
sensitive area planned for the balloon flight instrument. The detector is currently being tested for its cryogenic, 
spectroscopy and imaging performance. Here we will discuss the design philosophy and describe in details this TPC 
system, including the cryogenics, inner readout structure, and electronics system. The initial experimental results 
obtained with this large prototype are also presented here. 

2. LXe-TPC 10 LITER PROTOTYPE 


2.1. Xenon Purification System 


The purity of the LXe is an essential requirement for the operation of the proposed LXe-TPC. Fig. 2 shows 
the schematic of the purification and gas handling system which has been built and used for filling the large volume 
TPC. It uses commercially available purifiers and only UHV components. The purified gas is stored in two stainless 
steel high pressure cylinders which are also used to recover the gas from the TPC after each experiment. We have 
already demonstrated that a sequence of Oxisorb 12 , cold molecular sieves and high temperature getter is capable 
of purifying the xenon gas to the sub-ppb level. This purity level is needed to drift ionization electrons over large 
distances 7 . 

A gridded ionization chamber can be used to monitor the purity of the liquid xenon by detecting the signals 
of cosmic rays traversing the detection volume. The measurement is usually done at very low electric field strength 
to achieve long drift time, in order to get better sensitivity to the effect of election attachment by electronegative 
impurities. From a pulse shape analysis of the signals induced on the anode by the drifting electrons, information 
on the electron lifetime and thus on the total impurity concentration can be obtained. 

Fig. 3 is a typical cosmic ray signal registered in a 3.5 liter gridded ionization chamber used as purity monitor. 
A linearly rising signal with a long rise time is a direct indication of good liquid purity. An electron lifetime, longer 
than 1 msec in LXe, has been achieved. It indicates an impurity level of < 1 ppb O 2 equivalent. This purity level is 
more than sufficient for our proposed detector where the maximum drift time is not longer than 40 /xsec. 


2.2. Chamber and Cryostat 


The 10 liter LXe-TPC consists of an inner cylindrical vessel, in which we fill the ultra-pure LXe and mount the 
electrode structure, surrounded by an external vessel which is used as a cryostat. Both vessels are made of stainless 
steel. Fig. 4 shows the cross-section of the system. The volume between the two vessels is evacuated and maintained 
at a pressure of 10“ 5 ~ 10“ 6 Ton. To achieve better thermal insulation, a few layers of aluminized Mylar foil are 
used to surround the inner vessel. The inner vessel is constructed with standard UHV technique. All the construction 
materials are bakeable and capable of working at cryogenic temperature. This inner vessel is designed to withstand 
a pressure up to 5 atm. The normal operating pressure is around 1.5 atm. The inner vessel bottom is a 18-1/2" 
bolted flange, on which the whole electrode structure is supported. 

The cooling of the TPC is achieved by a controlled flow of LN 2 , through the copper coil of the condenser 
provided on top of the inner vessel. The coil has 1/4" diameter and ~ 15" effective length. The condenser is 
thermally insulated by a vacuum jacket up to the top, where the port for pumping and gas filling is located. During 
the operation of the system, the vapor pressure on top of the liquid is monitored by a pressure gauge and controlled 
by the LN 2 flow in the condenser. 

It has been tested that such a cooling system is very effective and efficient. It needs about 80 liters of LN 2 and 
4 hours to fill the 10 liter chamber, and about 5 liters of LN 2 every 3 hours to keep the operating condition. 
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Fig. 2 Schematic of the xenon gas purification and handling system 
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1 U TW °, V,Y, transparent windows on the inner bottom flange are optional for detecting the primary scintillation 
lg t wit U sensitive photomultipliers. Crystal quartz is chosen for its high transmission at the wavelength of 175 
nm of the xenon scintillation light. 

2.3. Electrode Structure 


For imaging the point-like charge blobs produced by MeV 7 -rays in LXe, we have proposed the non-destructive 
readout structure , schematically shown in Fig. 5. It consists of two orthogonal induction wire planes, followed by 
a plate anode for total energy measurement. The plate anode can be segmented into sections, each connected to 
its own charge sensitive amplifier, to reduce the input capacitance and hence the electronic noise. A shielding grid 
separates the sensing region from the drift region. A point-like charge produced within the sensitive volume and 
rifting to the wire planes has to pass through the shielding grid and the induction planes without loss, to be fully 
collected on the anode. This is achieved by choosing the proper electric field ratio before and after these planes. As 
soon as the ionization electrons are created and start moving in the drift region, they will induce a current pulse on 
the grid, which is a function of drifting distance and can be used to indicate the initial position where the charge 
was produced. Once passed through the grid, the electrons induce a positive current on the wires below, when they 
are approaching, and a negative current, when they are drifting away from the wires. Fig. 6 shows the current and 
charge signal waveforms on the grid, induction wire and anode produced by a point-like charge. While the induction 
signal is triangular in shape, the collection signal from the anode plate is step-like. The signals induced on the wires 
provide the X — Y coordinate information. The measured drift time, referred to a time zero, together with the 
known drift velocity, provides the ^-coordinate information. For 7 -ray events with multiple Compton interactions 
the height of each step in the resulting anode signal gives the energy deposited in each point, and the time interval 
between each step gives the relative distance of the interaction points, along the drift direction. The sum of all the 
step pulse heights is proportional to the total energy of the original 7 -ray, if totally absorbed. When combined with 
the additional X -Y information inferred from the induction signals, 3-D imaging and calorimetry is achieved for 


This readout structure also permits to improve the accuracy of the X — Y information, by using the point-like 
ionization feature produced by low energy 7 -rays in LXe. In fact, the amplitude of the signal induced by a point-like 
charge blob depends strongly on its initial lateral position with respect to the sensing wire. This position dependence 
allows to derive the spatial coordinate of each blob by weighting the signal amplitudes on neighboring wires. Thus 
the spatial resolution m the detector plane can be better than the wire spacing, and is determined by the accuracy 
of the amplitude measurement, i.e. the signal-to-noise ratio. Ultimately, the localization of a charge blob is limited 
by its physical extension and its diffusion during the drift. For a 1 MeV electron in LXe, the intrinsic size of the 
charge blob is on the order of 0.5 mm. The lateral diffusion, orthogonal to the drift direction, is given by y/DH 
where D = 65 cm s is the diffusion coefficient 13 , U is the drift time. For a maximum drift time of 40 nsec, as in’ 
the proposed LXe-TPC, the resulting value is about 0.5 mm 

This readout scheme has been implemented and tested with the 10 liter LXe-TPC prototype. Initial tests have 
been made with a 3 cm drift gap, to minimize the applied voltages on the grid and induction wires. 

The sensitive area of the anode is 19 cm x 19 cm. It is segmented into four identical 9 cm x 9 cm squares each 
coupled to its own preamplifier. The sensitive volume is thus only 1,200 cm*. Hie total liquid volume of 10 liters is 
however of the same order as the sensitive volume of the proposed fli ght instrument. 

The shielding grid is a mesh made of 100 jon diameter stainless steel wires with a pitch of 2 /im, orthogonally 
stretdied on the same ceramic bars on which the induction wires are stretched, to allow for good alignment and 
therefore minimum electron trapping. The sense wires are also made of 100 /im diameter stainless steel wires 
with 4 mm pitch for both X and Y readout. The gap between the different layers is also 4 mm. The mesh iri 
this arrangement has a shielding inefficiency of about 3%, according to the Bunemann’s formula 14 , modified for 2- 
dunenffloiwl grid. The fields above (£i) and below (Ej) the shielding grid are chosen to permit complete transparency 
to the diming electrons. Here, a field ratio E^/E\ larger or equal to two is used for 100% transmission. 

The assembled electrodes structure is mounted on the bottom flange. HV is separately provided to the cathode 


6 




Fig. 6 Signal waveforms on grid, induction wire and anode plate: a) initial charge position, b) grid wire, c) 
induction wire, d) anode plate 


and the grid. Commercially available multi-pin feedthroughs (10 pins each, rated for 3.3 kV), also mounted on the 
bottom flange, are used for the induction signals. 

2.4. Flront-End Electronics 


The requirements for the front-end electronics include low noise, low power consumption, and compactness. The 
“““J 1 ™ P^uced by 1 MeV radiation in LXe corresponds to ~ 60,000 e~. At a drift field of 1 kV/cm, only 
~ 80% of these electrons are collected. In order to view events with multiple Compton interactions, the X - y wires 
6houid be able to sense point-like charge blobs with a m i nimum energy deposition of 100 keV, or 5 000 e~ in LXe 
In addition the signal on the induction wire is less than the total charge and depends on the lateral distance from 
he charge blob to the wire. To obtain a reasonable S/N ratio, the noise from the preamplifier, with the HV applied 
has to be reduced to a level of ~ 1, 000 e~ rms or better. ’ 


The total number of sense wires for the multiwire readout structure under test is 96. Each sense wire is connected 
to a low noise charge sensitive hybrid FET preamplifier via a high voltage capacitor, to decouple the input from 
the apphed voltage on the wire. The signal decay time constant is 220 /nee for the induction wires and 1 msec for 
the anodes. All the preamplifiers are enclosed in the cryostat which provides a perfect shielding as a Faraday cage 
The power consumption of the preamplifiers should be as low as possible, given the power limitations imposed by a 
balloon flight experiment. The low power dissipation is also important to avoid the bubbling of the liquid, and hence 
the low frequency noise, due to heat transfer via thermal conduction and thermal radiation. 

With a bus of ±6V the present hybrids dissipate 50 mW/channel. The noise level, with the preamplifier 
connected to the wire is about 400 e~ rms, and slightly larger with the HV on the electrodes. 

The preamplifier outputs are brought outride the cryostat through multipin connectors mounted on the wall 
of the cryostat, and given to the second stage amplifiers via shielded twisted-pair cables for filtering and further 

amplification needed by the waveform digitisers. The shielded twisted pair cable can effectively avoid noise pick-up 
and cross-talk between channels. r r 
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2.5. Signal Triggering and Data Acquisition 


To test the imaging performance of the 10 liter LXe-TPC, 7 -ray data have been taken with the detector in a 
self- triggering mode, using the charge collection information recorded on the anode with external calibrated 7 -ray 
sources. In order to insure that the incoming 7 -ray energy is totally contained within the sensitive volume, we select 
an energy window around the full energy peak of the source spectrum. As shown in Fig. 7, the ionization signal from 
the anode is digitized, after proper amplification, as well as used for event triggering. To be used for triggering, the 
anode signal is passed through a shaping amplifier before being fed into a single channel analyzer which generates 
a logic gate signal whenever the pulse height corresponds to the selected full-energy peak. This logic signal, after 
proper delay, is used as a common stop for the CAMAC based LeCroy 2262 FADC system. 
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Fig. 7 Block diagram of front-end electronics and data acquisition system 


A distinguished feature of LeCroy 2262 FADC is its poet-triggering capability. This is needed for the present 
set-up, once the anode signal is detected later in time than the induction signal on the wires. The FADC has a 10 
bit resolution and a sampling period of 200 ns, which corresponds to 320 sampling points over a period of 64 pace. 

The digitised anode and induction wires waveforms for each 7-ray event are transferred to a ftVAX for online 
display and offline analysis. 

The TPC spectral response is directly measured from the anode signal, using a multichannel analyser (MCA). 
The anode signal is first amplified and shaped with a research amplifier or a gated integrator 11 . 

3. EXPERIMENTAL RESULTS AND DISCUSSIONS 

We have carried oat repeated experiments with the 10 liter LXe-TPC prototype and have achieved the primary 
d e si gn objectives: (1) verify that the 10 liters of liquid xenon can be purified and kept dean for good charge collection 
over a period of several days; (2) demonstrate simultaneous spectroscopy and imaging of MeV y-rays, with good 
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energy resolution and submillimeter spatial resolution; (3) use the measured energy and spatial information for event 
reconstruction by Compton kinematics and identification of the first Compton interaction point. 

In the following, we summarize the experimental results achieved so far. For the imaging performance, only 20 
of the induction wires were used, due to the available number of second stage amplifiers and digitizers. 

3-1. System Stability 


The long term stability of the LXe-TPC performance has been monitored for a running time of four days. The 
full-energy peak position of the 511 keV 7 -rays from a 22 Na source was used as an indicator of the charge stability 
as a function of time. Results are shown in Fig. 8 , indicating no appreciable decrease in charge collection over the 
entire period of the run. The collected charge is consistent with the result expected with good liquid purity. The 
fluctuation around the straight line fit is attributed to the error in peak position reading from different measurements 
and to the stability of the electronics and the high voltage power supply. A better control of the operating conditions 
should reduce this fluctuation. 



Time (hour) 


Rg.8 Charge collection yield from 5 1 1 keV r*ay as a function of time* obtained with 10 liter LXe-TPC 
at a drift field of 1 kV/cm 


3.2. Spectroscopy Performance 


ener ® r re 8 P° n » e WM measured with collimated m OB and M Na 7 -ray sources as a function of applied electric 
field in the range of 0.5 ~ 1.75 kV/em. typical pulse height spectra are shown in Fig. 9 for a drift field of 1 kV/cm. 
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The spectra were obtained with a MCA after a 6 
pulse distribution, indicating an electronics noise 


psec gated integrator shaping time. The rightmost peak is the 
contribution of about 400 e~ rms, or 20 keV FWHM. 


test 



Fig. 9 (a ) 137 Cs and (b) spectra obtamed with 10 liter LXe-WCatadrift field of 1 kV/an 
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ener gy resolution obtained with the 10 liter LXe-TPC is compared with that previously measured 
with simple gridded ionization chambers of much smaller volume and only a few millimeters of drift gap 9 . The open 
squares are the 570 keV FWHM energy resolution data, normalized to 662 keV by the E~ l/2 law. The closed circles 
fu/bm ata k° m ^ IC ^ lter LXe-TPC irradiated with 137 Cs. The noise-subtracted energy resolution is 12.1% 
FWHM at 662 keV, comparable to that in the much smaller chamber, at the same drift field of 1 kV/cm. This 
imp les that the liquid purity is sufficient for long drift and the signal processing with the gated integrator prevents 
the resolution from degradation due to the ballistic deficit as discussed in ref. 11. 


These results are very encouraging. Better energy resolution can be achieved at higher drift fields, by reducing 
* n electron - ion recombination. For example, at 4 kV/cm we have measured a resolution of 
7% FWHM for 662 keV 7 -rays with a 3.5 liter gridded chamber and 5 cm drift gap 16 . However, in the TPC mode of 
operation, the detection of small induced signals in presense of HV, makes this approach difficult. A more practical 
approach is provided by the photoionization effect in LXe doped with appropriate organic molecules, to convert the 
scintillation light to charge. The increased charge collection and reduced fluctuation in recombination results in 
improved energy resolution 15 . We started our studies in this direction with LXe doped with TMA (trimethylamine) 
wffich is believed to be the most efficient photoconverter with LXe. We have tested the spectroscopy response of 
LXe-TMA with a gridded ionization chamber of 1 liter volume. 


Fig. 11 shows the system used for the doping experiments. Xenon gas is purified by an Oxisorb followed by a 
spark purification chamber . TMA is purified separately with vacuum distillation and trapping with cold molecular 
sieves The test chamber includes a UV sapphire window to couple a PMT to the sensitive volume and study the 
quenching of the light by the dopant. We have observed encouraging spectral improvement with even a small 
concentration of TMA (7.5 ppm). Fig. 12 shows the spectra of 207 Bi obtained in pure xenon and doped xenon at 

Tin^v d °! 5 kV / cm - The better separation of the two higher energy lines (976 keV conversion electron 
and 1048 keV 7 -ray) can be clearly seen in the spectrum with doped xenon. In Fig. 13, the lower curve is the field 
dependence of the energy resolution in LXe with 7.5 ppm TMA, while the upper curve is obtained with pure LXe 
lo test the stability of the detector’s response with time, we operated the chamber for more than 7 hours No 
appreciable change m charge collection was observed (see Fig. 14). We will continue these experiments with different 
concentrations of TMA to find the optimum in energy resolution and system stability. 



Fig. 11 Test chamber and purification system for studies with TMA doped LXe 
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FWHM Energy Resolution (%) 



Channel Number Channel Number 


Hg. 12 Comparison of 207 Bi spectra obtained with (a) IMA doped LXe (b) pure LXe at tbe same drift 
field of SkV/cm 



Hg. 13 Field dependence of FWHM energy re soia tiop 
for 570 kcV y-ay in (o) pure LXe and (•) LXe doped 
with 7.5 ppm TMA 


Fig. 14 Col l ected charge as a fimctiou of time for 
570 loeV y-cay interactfon inside LXe doped 
with 7.5 ppm TMA 
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3.3. Imaging Performance 


Due to the limited number of readout electronic channels, we used only part of the readout structure (10 X- 
wires, 10 Y- wires, and the corresponding section of the anode plate located below these sense wires, corresponding 
to a sensitive area of 4 cm x 4 cm) to obtain the spectrum and imaging of a collimated 7-ray beam from 137 Cs. The 
collimator is made with two lead blocks with a thickness of 5 cm and a hole of 2 mm diameter. 



Fig, 15 2~D distribution of -2000 events of 662 Ice V Y-rays from collimated 137 Cs source detected with 
the 10 liter LXe-TPC. The diameter of the collimator is 2 mm 


By weighting the signals induced on two neighboring sense wires by a charge blob produced by a 7 -ray one 
can locate the interaction point coordinates in the X - Y plane with a precision much better than the 4 nun’ wire 
spacing, triggering on 662 lteV 7 -mys from the full-energy peak, the collimated 2 mm diameter beam was imaged by 
this readout structure. Fig. 15 shows the 2~D image of this vertically incident beam. The beam width is measured 
as = 1.6 mm, including the possible spread due to the solid angle effect of the collimation, the Q f the 

beam, as well as the spread due to Compton scattering. The centroid finding algorithm used for this analysis needs 

to be optimized. More data have been accumulated and we expect that a cleaner image with less background will 
be obtained. 

The preliminary imaging results of Fig. 15 together with the good energy resolution shown in Fig. 10, obtained 
simultaneously, represents the realization of the second design objective of this LXe-TPC development. 


3.4. Compton Ki n e ma t ics and Event Reconstruction 


It is another design objective and a unique capability of the LXe-TPC to identify the first Compton interaction 
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point and to reconstruct the incident 7-ray direction by means of Compton kinematics. The former is required 
by the coded aperture imaging and the latter is needed to identify and reject background. The following example 
demonstrates how the measured energy and spatial information for a Compton scattering event is used to reconstruct 
its topology. Combining the induction signals on the X - Y plane with the collection signal on the anode the (X 
K , Z) coordinates as well as the energy deposition E at the scattering and absorption points can be obtained With 
this information, Compton kinematics can be used to find if the incident 7 -ray direction is within the FOV Fig 16 
is the on-line display of a Co 7 -ray event, showing the digitized signals from the 10 Jf-wire, 10 Y-wires and the 
ano e. The source was placed 10 cm above the 4 cm x 4 cm sensitive area of the detector, to approximate a celestial 
7-ray point source. The event was recorded with the trigger from the 1.33 MeV full energy peak of 60 Co. Channel 
1 ~ 10 correspond to the A"-wires, 11 ~ 19 to the Y-wires and channel 20 corresponds to the anode waveform The 
zoom of channel 20 clearly shows a 2-step pulse shape. The amplitude of the steps is proportional to the energy 

deposited in the two interaction points, their sum corresponding to the total energy of 1.33 MeV. The time interval 
between the two steps gives the vertical separation Z of the two interaction points, given the known drift velocity at 1 
kV/cm. Each step has its corresponding induction pulses on two neighboring .Y-wires and two Y-wires, respectively. 


We define the coordinates of channel 1 and channel 11 as the origin of the X — Y plane (X = 0 and 
respectively) and the first charge blob in time is at Z = 0. From the analysis of the waveforms in Figure 
two charge blobs (A and B) have coordinates and energy depositions as follows: 


y = 0 , 

14, the 



X (mm) 

Y (mm) 

Z (mm) 

E (MeV) 

Blob A 

20.1 

23.2 

0.0 

0.60 

Blob B 

28.0 

28.0 

8.0 

0.73 


1 . 2 wa f inferred from the drift time interval between two steps in anode signal, X and Y coordinates were 

obtained by weighting the amplitudes of the induced signals, and E is the energy deposition at the interaction site 

ova ‘° i * 1 pote Tke " ^ ; 

.. f g A . 17shows ‘ he reconstructed event in 3-D space. The most probable incident direction is vertically down along 
the Z direction. From geometry, the angle 0, is 49°, obtained from tan*, = f/Z, where / is the lateral separation 

bteT^Lnpdnt ^ 4186 ^ 8Catterin * an « ,e * 6 <' be spared to*, to directly identify the first 

Eq 

= ii 3 Me , V ' ***** of the incoming 7 -ray. If the scattered photon has energy E 7 = 0.60 
MeV, then 0 C - 58 , and Compton kmematics is not satisfied. If the scattered photon has energy Eh = 0 73 MeV 

the first Compton interaction point and the event as a case of forward scatterin g 

P° mpt ™ e 5 ve * “T***?*? procedure can be generally done with a reconstruction algorithm based on 
Compton kmematics . Following the formalism of ref. 5 , we define the function f as: 

f=Y,W-W!) 2 

«= o 

where is the energy of the ith scattered photon calculated from the Compton formula, while W! is the enersv of 

toSl 8CatterC<1 P f t ° n J; Ut COnServation on] y- To test the validity of the assumption That 

the total energy is contained, as well as if the assumed initial direction of the 7-ray is kinematically possible this 
function is m i nimiz ed for all possible combination of points. ’ 
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For the two-site event in Fig. 17, the function / has only one term 


/ = W - Wif = - (£, + g,)] : 


-hid, gives the smaller / value. / = 0 23 if 2 -S.fl.tT "« -^tilied by the energy assignment 
interaction point. ‘ ~ ° 73 M ' V ' but ! = 0 01 lf E ‘ = MO MeV. Blob A is the first 

for ^e^^Wve^udumT^An Co* 11 ^ 7 

reconstruction of the direction of the incoming -y-rav and a LXe-TPP . ’ ^“P 4011 ^nematics enables the 

in the original paper 2 . However, the resultingtlgulw resolution for th^M v** " aC ^ m P ton teIe8co P e > “ discussed 
between two successive 7-ray interactions in the dense LXe While not* r Iuiuted by the 8maU separation 

angular resolution is however good enough for e^ent^ l~ «■ ^ 

mds £ '£xzi££i TJZSSUift: ™i' orf “i' rf *r; ofl “ lf " gl ' s ‘'“ d ' 1 ''^ 

deposition at different interaction points using the Compton kiLmafuT^ n ftora * he me “ u «d energy 

cone onto the celestial sphere is tiled the event ci^r ^e^rt'hTdir 0 ^ 0 ™' ?! ProjKtio ° ° f ,bis 
the scattered y-ray, and the event cirde radius is identical to the derived r t“ determined by the direction of 
sonrce is then located as the intention pom. of aff dm individud event P ° int 

r * ” n * -*■* — »•— *> - 

we covered the sensitive area by a lead blik^rith a 1 cm “*}*’, an l thus approach a parallel flux situation, 

two interaction points were selected. The Compton imaging capacity rf^rxScf^b 110 ^ ****** Wlth 
event circle outside the FOV can be identified and w 1 , LX e- Tp C is obvious. Any event with 

reconstruction based on Compton kinematics does not require the ahaS 0 ^ ** WC h&Ve shown > the event 

provided by the light trigger. However, the S*Tte Z w ^ W * *• **0 

energy background events by fiducial volume cuts 5 Without livht t *” mtera f tlo “ P° lnt “ necessary to reject low 
this information can still be obtained from the prompt grid sWl ““ ° f *, d ? ed LXe detector - 

signal as recorded for a 7-ray event in our TPC The «t*fw^ ' eth 19 J* ows for «ample the grid and anode 
the interaction point. our TPC. The starting t.me of the grid signal directly gives the Z location of 


1 . CONCLUSIONS AND P ER-SPF.riTTVi?. 


A new type of 7-ray telescope has been proposed to imare the MeV ~ „„„ a , , . 

sources, with high efficiency, good energy riolluhm ^ 7 T* flUX &0m galactlc “ d extragalactic 

lx - tpc - » *— -*«- sszJi: 

” d ve*> puii^ff lmStottuh^klh^ SST l “ b<! “ W “ 

“ u “ ^ « - « - * iv/cm, to JSftrs.ti: 

^^Impmvumtu. of thu tuutg, purlurm^of Ufi, mmgi,gd«u«m mil b« «hiu,d by mumfTMA dopud 
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Fig. 18 Reconstruction of the "event circles" for 22 “Na y-ray events from a point source 
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Kg. 19 Grid and Anode signal for a ^Na y-ray interaction in the LXe-TPC. The time scale is 5 ps/div. 
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.edJ^r^rir' UliOT "" be ^ “* »f K eiectromm au d by 

The ultimate spatial resolution limit, both for the Z and the X _ v t 
fin,te “‘ en5i0n of th ' b '° b i» which is on the order ot 0.5 mm fof™ MeV eLctom °" "" ‘ h ' 
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